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We present a time-dependent calculation of the velocity-dependent forces which act on an excited 
atomic dipole in relative motion with respect to ground state atoms of a different kind. Both, its 
interaction with a single atom and with a dilute atomic plate are evaluated. In either case, the total 
force consists of a conservative van der Waals component and a non-conservative Rontgen compo¬ 
nent. On physical grounds, the former corresponds to the velocity-dependent recoil experienced by 
the excited atom in the processes of absorption and emission of the photons that it exchanges with 
the ground state atoms on a periodic basis. The latter corresponds to the time-variation of the 
Rontgen momentum, which is also mediated by the periodic exchange of quasi-resonant photons. 
We find that, at leading order, all these interactions are linear in the velocity. In the non-retarded 
regime the van der Waals force dominates, being antiparallel to the velocity. On the contrary, in the 
retarded regime the velocity-dependent forces oscillate in space, van der Waals and Rontgen forces 
are of the same order in the atom-atom interaction, and the Rontgen component dominates in the 
atom-surface interaction. 


I. INTRODUCTION 

Frictional forces mediated by the electromagnetic 
(EM) vacuum field between neutral atoms and dielec¬ 
tric objects in relative motion has been the subject a 
number of recent works [IHI3- As already appreciated 
by Einstein in his celebrated paper of 1917 [il|, it is 
the Doppler shift on the photons which interact with a 
neutral atom in motion with respect to a thermal bath 
that gives rise to an effective friction which drives the 
atom to thermal equilibrium. Since then, a number of 
authors have discussed the possibility of an analogous ef¬ 
fect at zero temperature fll-llOl fia - ll^ . In this case, it is 
the exchange of vacuum photons between the fluctuating 
currents of neutral objects moving at relative velocities 
that gives rise to a net transfer of momentum between 
them in the direction parallel to the velocity. In this re¬ 
spect, the explanatory work by Pendry [ij has motivated 
the investigation of the so-called quantum friction in dif¬ 
ferent scenarios. Several approaches and approximations 
have been considered since then by a number of authors, 
yielding apparent contradictory results with regard to the 
dependence of the friction force on the velocity. Nonethe¬ 
less, recent works agree in the cubic scaling of the friction 
force with velocity on a ground state atom at zero tem¬ 
perature in motion with respect to a metallic slab imsi- 
For the case of the interaction between two macroscopic 
objects the approaches are semiclassical and based on lin¬ 
ear response theory and the fluctuation dissipation theo¬ 
rem ii. Concerning the interaction between an atom 
and a dielectric object, the authors of Ref. [1,11, nil have 
treated the atom quantum-mechanically, while the elec¬ 
tric response of the system macroscopic object-EM field 
has been considered semiclassical. On the contrary, the 
authors of Refs. HiiSlil have applied time-dependent 
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perturbation theory within a Hamiltonian approach in 
which the interaction between the macroscopic objects 
and the EM field is effective and semiclassical. Lastly, the 
interaction between two harmonic oscillators (intended to 
describe two atoms) has been considered by the authors 
of Ref. [1, B, [III within a Hamiltonian approach in which 
the EM interaction between them reduces to the electro¬ 
static (i.e., non-quantum) one. Finally, all the aforemen¬ 
tioned works, except for that by Scheel and Buhmann 
[1|, deal with bodies in their (internal) ground states, 
and they all appeal to some effective mechanism of dissi¬ 
pation to explain quantum friction. 




FIG. 1: (a) Sketch of the interaction of atom A with an atom 
B. Atom A moves at constant velocity v while atom B re¬ 
mains at rest at a distance R which varies in time, (b) Sketch 
of the interaction of an atom A with a thin plate made of a 
random distribution of atoms of kind B, with numerical sur¬ 
face density a. The atom A flies with velocity v parallel to 
the plate at a height d. The wavy lines represent the pairs of 
photons exchanged by the atoms. Emission and absorption 
processes take place at different locations of atom A. 

In this article we study the velocity-dependent EM 
force between two neutral bodies, one of which is ini¬ 
tially prepared in an excited state. In the first place 













2 


we study the interaction between two atoms of different 
kinds, A and i?, in relative motion, one of which, say 
atom A, is initially prepared in an excited state at zero 
temperature-see FigUJa). This is probably the simplest 
realistic scenario where a Hamiltonian treatment can be 
carried out all the way to the end. Important is the fact 
that dissipation is purely radiative here, which allows us 
to trace it along the calculation. In order to simplify the 
calculation we adopt the quasiresonant approximation 
outlined in Ref. M- That is, we assume two-level atoms 
for which the detuning between their resonant frequen¬ 
cies, Aab = — is such that Ta.b < Aab “C uja,b, 

where Tyi ^ are the linewidths of the corresponding tran¬ 
sitions of atoms A and B respectively. This condition 
is easily met in pairs of alkali atoms as well as, in an 
approximate manner, in pairs of circular Rydberg atoms 
[T^ . It allows us to truncate our results at leading order 
in Aa.b/^ab- In addition, we assume that the obser¬ 
vation time T satisfies 2'k/\Aab\ < T <C 2'kITa,b^ such 
that atom A remains excited until the obervation time. 
We will extend our approach to the calculation of the 
velocity-dependent force on an excited atom of kind A 
flying at constant velocity parallel to a thin plate made 
of a random distribution of atoms of kind B in their 
ground state -see FiglTJb). 

The article is organized as follows. In Secln]we explain 
the fundamentals of our approach. In Sec lIIIl we compute 
the velocity-dependent van der Waals force. The Rbntgen 
force is evaluated in Sec lIVI Conclusions are summarized 
in Sec El 


II. FUNDAMENTALS OF THE CALCULATION 


We aim at computing the velocity-dependent force on 
the atom A as a result of its EM interaction and its rel¬ 
ative motion with respect to atom B, in the electric- 
dipole and quasiresonant approximations. Our approach 
is based on that of Ref. |1^ . which is appropriately ex¬ 
tended to account for dynamical effects. Basically, the 
extention consists of, in the first place, promoting the 
center of mass (CM) degrees of freedom (d.o.f.) of the 
atoms to quantum operators, and secondarily, adding 
the contribution of the Rbntgen momentum both to the 
forc e op erator and to the interaction potential. Following 
Ref.|I 8 j|. in the electric-dipole approximation the canon¬ 
ical conjugate momenta of the atomic CM position vec¬ 
tors are Pa,b = Qa,b — dn,B ■ Qa.b 

being the kinetic momenta of the CM of each atom. The 
interaction potential between the free EM field and the 
atomic dipoles is IF = Wa +W§ + Wa + W§, where 


W^B 

wIb 


—dn.B ■ E(Ra,_b), 

Pa,b ■ [dA,B X B(Ra,_b)] 


+ 


B(Rn.B)] • Pa.b l2mA,B 


are the electric field and Rbntgen coupling potentials 
respectively which act upon the free Hamiltonians. 
These are, the free Hamiltonian of the atomic inter¬ 
nal d.o.f., Ha + Hb = hujA\A+){A+\ + hujB\B+){B+\, 
the kinetic energy of the atomic CMs, ]Ca,b = 

c^Q\ g + c^m\ g, and the Hamiltonian of free photons. 

Hem = e e + 1 / 2 )- I^ til® above express- 

sions Aa^b denote the electric dipole moment operators 
of the atoms A and B respectively, mA,B are their re¬ 
spective atomic masses, and E, B are the electric and 
magnetic field operators which can be expressed as series 
of annihilation and creation field operators of photons of 
momentum fik, respectively, 


E(Rn,B) = Y.^i~A'^AB)+W\riA,B) 

k 


E l hck , 


2Vei 


0 


ik-RA.B _ g —ik-RA.Bj 


B(Ra.b) - ^Bi^'^(RA.B) + B^+)(RA,B) 


= *E 


k,e 


2cfcVeo 


k X [eak,£e 


ik-RA.B _ c*£.f „-*k.RA,. 


where V is a quantisation volume, ui = ck is the photon 
frequency, and the operators aj^ ^ and Ok.e are the cre¬ 
ation and annihilation operators of photons of momen¬ 
tum Kk and polarization e respectively. 

As in Ref. m , we use the wave function formalism 
in the Schrbdinger picture [l9j | to compute the veloc¬ 
ity dependent-force on an excited atom A which moves 
at certain velocity v with respect to (w.r.t.) an atom 
B in its ground state. That is, denoting by U(r) the 
time-evolution operator at the time of observation T, 

ILl(T') = H exp —ih ^ Jp At [Ha Hb Hem + + 

ICb + W^ + W§ + W^ + W§]Y and by |4'(0)) the wave 
function of the two-atom system right at the moment at 
which atom A gets excited, the total force acting upon 
atom A at r > 0 reads, at 0[{Wa + IF 4 )^], 


(Fa(T)) = driQAiT)) = -zIi5T(^(0)|Ul(r)VR^U(T)|^(0)) 

+ aT(«'(0)|ut(T)dA X B(R^)u(r)|4'(0)) 

-VR,([IFf + W^]{T))/2 + aT([d^ x B(Rn)](T)) 

~ ivR^(dA-[E(R^) + m^i(Qn) xB(R^)](T)) 

+ dT([dAxB(R^)](r)). ( 1 ) 


In this equation a term quadratic in B(Ra) has been 
discarded in the last line as it yields a negligible con¬ 
tribution under quasiresonant conditions and we as¬ 
sume that quantum fluctuations over the classical tra¬ 
jectories of the atoms are negligible, both in the posi¬ 
tions and in the kinetic momenta of their CMs. That 
is, denoting by Ua:(T) = exp {—+ /Cs]} the 
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free-evolution operator of the CM d.o.f., the (approxi¬ 
mately) constant values of the atomic kinetic momenta 
read (Qa.b) = (^'(0)|U^(t)QA,BUK:(i)|^'(0)) V t. Next, 
neglecting quantum fluctuations over the expectation 
values of the CM position and momentum vectors, 
we can write (4'(0)|l[jj(-(t) J'(Ra,b, Qyi,B)UAc(i)|^(0)) ^ 
b)) for any functional J- of the opera¬ 
tors and Qa,b- In particular, the functional gra¬ 

dient Vr^ is replaced with Vr,, in Eq.(IT]), with Rt = 
([Ra — RR](r)), and the constant velocity of atom A re¬ 
lates to (Q^) through the classical relativistic expression 

V=^(»(o)nipi)R,i]^(i)i»(o))= 

Similar to the calculation carried out in Ref. [l6l| for the 
van der Waals (vdW) interaction, we will apply standard 
time-dependent perturbation theory at order 4 in kC to 
compute the force of Eq. m under quasiresonant condi¬ 


tions. The only difference here is the additional coupling 
to the Rdntgen momentum. We will address first the 
velocity-dependent force in the absence of the Rontgen 
momentum, which we will refer to as velocity-dependent 
vdW force. Later, we will compute the terms of the force 
which depend on the Rontgen momentum, which we will 
refer to as Rdntgen force for brevity. 


III. VELOCITY-DEPENDENT VDW FORCE 

It was found in Ref, [l^ . using standard time- 
dependent perturbation theory at order 4 in the potential 
+ Wg that the energy of interaction between an ex¬ 
cited atom A and a ground state atom B at time T after 
the excitation of atom A is given by 


J 






vPdk /■°° Vk'^dk' 

^—rrr / ^—tt;— / dil di2 

(27r)3 (27r)3 Jq 


oo v / —oo 
T rt rt' 


z(^(O)|Uo(-r)|4-(O))0(r - 2R/c) 


X f dt [ dt' [ dt"(^( 0 )|d^-E^ 7 ^(R^)Uo(r-t)ds-E[^)(RB)Uo(t-t')dB-E^"^(RB) 

Jo Jo Jo 

X Uo(t' - OdA ■ E'^+^(R^)Uo(t")|^( 0 ))] +[k^ k']\ 

I- 


( 2 ) 


where R the interatomic displacement vector, R = R^ — 
Rr and lUo(t) is the unperturbed time-evolution propa¬ 
gator in the absence of W. Eq.([2]) corresponds almost 
completely to the diagram of Figl^Ka) alone, where the 
photonic loop is made of two doubly-resonant photons. 
In addition, the diagrams (5) — (/) of hgure 1 in Ref.fl^ 
provide the Heaviside function 0(T — 2R/c) which guar¬ 
antees causality. 


For both atoms at rest (which we denote by the super¬ 
script 0) their position vectors Ra and Rr can be con¬ 
sidered classical and constant vectors, and Uo{t) reads 
U[)(f) = exp [-ih~^{HA + Hb + HEM)t]. Finally, |4'(0)) 
at zero velocity is made of the tensor product of the EM 
vacuum state (|0-^)) and the internal states of the atoms, 
the excited one for atom A and the fundamental one for 
atom R, |4'°(0)) = |H+) (g) |R_) ® |0.y). Substituting next 
these formulas in Eq.(I2]), taking the limit to the contin¬ 
uum by replacing the sums over momenta by integrals, 
summing over polarizations and integrating in orienta¬ 


tions we arrive at [1 

^3 


(Wf(T))° dkk^lm{G^^ik'R)} 


TTCoe 

^+oo 


dfc' /c'Hm{GP«(fci?)}0(r - 2R/c) 


(3) 


/o 


T i-t i-t’ 

dt / dt' / dt" 
Jo Jo 




-p (w o w')! 


where af is the fine structure constant, = 

(H_|d^|H+), = (R_|dR|i?+), and the tensor G{kR) 

is the dyadic Green’s function of the electric field induced 
by an electric dipole of frequency ck in free space, 

k 

G{kR) = — - [a/kR -t ij3/[kRf - l3/{kR)% (4) 

47r 

where the tensors a and (3 read a = I - RR/R^ /3 = 
I - 3RR/R2. 

Let us consider next that the atoms move away from 
each other at relative velocity v. For simplicity, we con¬ 
sider atom B at rest in our reference frame. We will 
denote the associated observables by the superscript v. 
At leading order, the dominant contribution to the in¬ 
teraction still comes from Eq.(I2|), but now allowing for 
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FIG. 2: (a) Diagrammatic representation of Eq.([2|) for 

(W’J'(r)). Time runs along the vertical as indicated by the 
arrows. R is intended, generally, as a quantum operator. 
Only for both atoms at rest, R can be treated in this picture 
as a classical variable, (b) Diagrammatic representation of 
Eq.® for {Wa{T))'' up to order (w/c)^. The interatomic 
displacement is here a classical variable which depends on 
time. v(r — t”) is the displacement of atom A during the lag 
between the emission of the first photon and the absorbtion 
of the second photon. 


the constant motion of the CM of each atom. In order 
to account for this motion, the position vectors of the 
CMs, R. 4 ,_b, and their kinetic momenta, Qa.b, are pro¬ 
moted to quantum operators; the initial state is now the 
tensor product of the EM vacuum, the initial internal 
atomic states and the states of the CMs of each atom; 
and the time-evolution propagator contains the kinetic 
energy of the CMs defined earlier, Uo(t) = UQ(t) ®UAc:(t). 
When the CM d.o.f. are considered quantum operators 
in Eq.Q, the resultant equation reads 

{wf{T)y c, 

J dtdt'dt"dkdk'[{'i'{0)\C{A±,B±,r, k, k', t, t", T) 

X U^(t' - + [A o k'W ( 5 ) 


where we have isolated explicitly the factors which de¬ 
pend on the CM operators alone, with C being a func¬ 
tion of the photonic and internal atomic d.o.f. as well 
as of the time and frequency variables. In the absence 
of the Rontgen momentum, the following relation holds 
in application of the canonical commutaion relations, 


g-ik.RA,Bg-iK(QA,B)t/?!.g*k.Ryi,B _ ^-iK{QA,B+Kk.)t/h 


X exp 


mA,B 


Qa,bI‘^ 


{mA,BcY 


( 6 ) 


In the last line we have expanded K, up to terms of order 
5 in Qa,bIbia,bc. Applying the above relation to Eq.® 


we end up with the expression, 


{Wf{T)y ~ J dtdt'dt"dkdk' 

(-t')C(A±, B±, 7 ; k, k', t, t', t", T) 


X exp 


X exp 






niA ^ 2{mAcY 8(mAc)^ 


-VO ■ 

2{mBcY ' 8 ( 77130 )^^ 


A QS ,,Nr- Qb 


X g®k''(R-A —Rslg — *k.(R. 4 — Rb) 


X exp 


Q-4 /,/ ^A 


3Q\ 


TTiA^ 2{mAcY ' 8{mAcY‘ 


xU 7 c(t')l^( 0 ))} + {fcofc'}t, 


(7) 


where terms of the order of hk/Q a.b times less have been 
ignored in the exponentials [^. We observe that the 
factors which depend on the CM d.o.f. are evaluated on 
the state of the CM at time t', UK:(tOI'I'(0))- Neglecting 
quantum fluctuations in the CM d.o.f. and making use 
of the approximation 


(Qa) L _ {Qa? {Qa)^ ' 
tua 2{mAcY 8(mAc)‘^_ 


( 8 ) 


we find 


(W^(T)y ~ j dtdtAf'dkdk' 

X [(vk(0)|C(A±, R±, 7 ; fc, fc', t, t', t", T)|vI/(0)) 

X g-»k.[<RA(T))-v(T-t")-{RB)]g*k'.[{RA(T)>-<RB>] (9) 

X -I- [fc o ky. 


This result shows that, up to terms of order 3 in v/c, 
the net effect of the application of Eqs. and m is 
the substitution of the operators Ra,b in Eq.Q by their 
corresponding classical vectors which evolve in time as 
{RbY)) = (Rb) constant, {RaY)) = {RA{T))-v{T-t), 
with 0 < t < T. The absence of Doppler shifts of order 
(v/cY is remarkable. 

The resultant diagrammatic representation of 
(Wj'(T))'' up to 0[{v/cY] is shown in Figl^Ib). Once 
integrated in orientations, the interaction energy of the 
moving atom A reads. 


{W^{T)y ~ f dA:' 


^+oo 


dfc 


TTCoe 

dt [ dt' 
0 Jo 


dt" 


k'ym{G^yk'R t)} 


X kym{GPyk\RT - v(T - t")|]} (10) 

X 0 [r(l -b vr/c) - 2Rt/c] {i 


X e 


— i{t — t')uiBp — i(t'—t")ui —it"uiA'\ 


) + [w -O' a;']'^|. 


where Rt was defined earlier as the interatomic displace¬ 
ment vector at the time of observation. In Ea. dTUl) the 
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argument of the first Green function is shifted w.r.t. the 
one in Eq. ([31) by an amount proportional to the displace¬ 
ment of atom A during the lag between the emission of 
the first photon and the absorbtion of the second photon 
by atom A, v(T — t"). Correspondingly, the argument 


of the Heaviside function is shifted by an amount vjiTjc, 
where vn is the component of v parallel to Rt- The shift 
in the argument of the Green function gives rise to two 
kinds of terms at 0{vr/c), 


G[/c|Rt- v(r-t")|] 


^ik(RT-VRT+VRt") - 


OCT 




Pt 


47r 


lk[RT - vr{T - t")] [k{RT - vrT - [k(RT - - t"))f i 


AkRx r 


47r 


-ikvR{T-t")^ Q;t iPr 

^kRr [kRrf 


Pt 


{kRrf 


+ 


vr(T — t") ar 


2iPT 


3Pt 


Rx 

I— 


^kRr [kRrY ikRrp 


( 11 ) 


where the tensors read ar = I — RtRt/Rt, Pt = 
I — 3RtRt/Rt- The terms within the first small square 
brackets of the second line of Eq. (EH) can be interpreted 
as Doppler shifts on the frequency of the photon emit¬ 
ted from atom A, whereas the terms within the second 
small square brackets are multiplied by the lag between 
the photons emitted and absorbed by atom A. This dis¬ 
tinction is however merely formal, since in both cases the 
origin of these terms lies in the Doppler shifts of the com¬ 
plex exponentials of Eq. m, restricted to terms linear in 
vrIc. 

The subtitution of Eq. EH in Eg. lfTOl) yields three dif¬ 
ferent contributions to (Wj'(T))'". We will denote the 
interaction associated to the first term of Eq. EH by 
and the interaction associated to the sec¬ 
ond term of Ea. lfTTl) by (Hbi(T))Xag' Lastly, the Heav¬ 
iside function contributes at leading order with a term 


(bb|)e = [21 ?t(1 - vn/c)/c])°U=o- In sum¬ 

mary, we obtain 


_ ‘^ijpq^ABVR \ gij ^pq d 2 


{wnm = 




_ k%RUP'^pP? + 2a^iP^^^) 


kgRPpa^aP^]^ sin — 


^CAijpq^ABTR 


RU 


kB 


X pl^P^ — k%R^a^P^ cos {2kART), 


( 12 ) 


which is of the order of vrAabRt/ c^ times the inter¬ 
action {Wa{T))° calculated in Ref.[l^, with Uijpq = 
fip /[{4:Tr€o)‘^hAAB]- As for the Doppler term we 

find 


{Wa{T))1,p = - ~k\Rl{P^P^T^ + 2 cP^P^t^) + k\RU^aPp^] cos (2^^^^) + ^^Ia^PI 

- k^RWipip] sin {2kART) - - klRUpf^pP? + 2 ab^) + k%R^a^a^^^] 


X cos 


{2kBRT + AabT) - ^^ksipt^P^P - klR^a^pi;?] sin {2kBRT + AabT), 


(13) 


which corresponds to the result obtained in Ref.jl^ for {Wa{T))° but for the replacement of and Aab with 
their Doppler shifted values, kA ^ kA = ^^(1 — vr/c), Aab = ^A — ^j^b —t Aab = ^A — wb(1 -|- vr/c) and 
Uijpq -A U^jpq = k-tk-q k-f k-p /[{^EtofhAAB]- Filially, the lag term reads 


Ka 


{WA{T))l^g = {vR/RTAAB)‘2lA,,pq{[2c-^AABRT cos {2kART) - sin {2kART)][3P^^P^TVR% + ‘^T^^PtPt + ‘^cP^Pt) 

Hrj 


+ k\cP^aP;C/Rl] - [2c-^AABRTAn(2kART) + cos {2kART)]{^k aP^P^VRt “ 3k\cP^P^^^/R^j,)], (14) 


where oscillating transient terms of frequency much 
higher than Aab have been omitted for simplicity. 


Next, we proceed to compute the velocity-dependent 
vdW force on atom A at order v/c, which is the 
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gradient of the difference between the electric interaction 
at velocity v ^ 0 and at zero velocity, 

{^AOiv)) = (T))Bop - «(T))° 

+ {W^{T))l + {Wf{T))l^^]. (15) 

Averaging in time Ea. lfT^ for T ^ averaging 

over dipole orientations and restricting ourselves to the 
near-field and far-field regimes, for which UaRt ^ 1 and 
UaRt ^ 1 respectively, we are left with 

UaRt « 1, 

(noT^)) ^ 4^^^«Nn(2fcAi?T) (16) 

— 2c~^AyiRi?^ cos ( 2 fc^i? 7 ’)], kARx ^ 


where U = /[{^Treo)^hAAs]- Interestingly, while 

in the near-field both Doppler and lag terms contribute 
at the same order to (F()(^^^), Doppler terms domi¬ 
nate in the far-field. Terms coming from {W^{T))q are 
negligible in both cases. In addition, while (F'^a^^v)) 
is antiparallel to vr in the near-field so that it can 
be referred to as a friction force, it oscillates in space 
in the far field. Physically, corresponds to 

the velocity-dependent recoil experienced by atom A 
during the processes of absorption and emission of the 
photons that it exchanges with atom B, which takes 
place at a rate |A^b|. ^'ao[v)) indeed propor¬ 

tional to the adimensional factor vxiuia,bIcAab-, and 
hence to the component of the decay rate of atom A 
which depends on the presence of atom B, = 

-EJrf'fcftk5^(l(^-i?-,7k..|U(T)|vh(0))|2) The 

expression within parenthesis is the probability of emis¬ 
sion of a single photon of momentum ftk and polarization 
e at a time T after the excitation of atom A. 

Finally, for the sake of completeness we compute the 
velocity-dependent vdW force on the excited atom A as it 
flies with velocity v parallel to a thin plate (11) parallel 
to the xy plane and made of a random distribution of 
ground state atoms of the kind B. The numerical atomic 
surface density is denoted by a, and the distance of atom 
A to the plate by d -see FigHKb). For simplicity, we 
adopt the pairwise summation approximation, according 
to which the total force over atom A is the sum of all the 
pair forces between atom A and each atom B. Equipped 
with our computation of the problem becomes 

simply geometrical, and its result is 


/ OO PCX) 

dx dy(F(i‘g^,)(x,y,d)) 

-OO j —OO 

U {kA + ks) , J , 

= A.,,, V, 

(r!S 3 I.))n - 

— 2c~^AABdCOS {2kAd)], ^ 1. (17) 


In contrast to the result of Ref. S’ is always 

a friction force for <C 1, while the direction of the 
force of Eq.(79) in Ref.Q depends on the sign of Aab- 
However, this essential discrepancy is due to the partic¬ 
ular form of the effective permittivity adopted by the 
authors to characterize the plate there -i.e., a Drude- 
Lorentz form, and not to an essential difference between 
their approach and ours. As a matter of fact, insert¬ 
ing the polarisability of atom B in the scattering Green’s 
function of Eq.(65) in Ref.Q and adopting the quasi¬ 
resonant approximation, a term is obtained of the same 
order as that in our Eg. lfTHl) in the non-retarded regime, 
with different numerical prefactors though. 


IV. RONTGEN FORCE 


In this section we consider the force induced by the 
coupling to the Rontgen momentum and the force due 
to the time variation of the Rontgen momentum itself. 
They correspond in Eq.(IT|), respectively, to conservative 
and non-conservative forces. 

As for the conservative force induced by the inclu¬ 
sion of the interaction in Eq. (El), it does not only 
yield an additional potential energy, (W ^), but also con¬ 
tributes to {W^) through the insertion of in the 
time-evolution operators. Since is already linear in 
(Qa), it provides terms of orders v/c, {v/c)^ and higher 
to the conservative force without the need of considering 
the quantum nature of the CM d.o.f. within the EM op¬ 
erators. Thus, an analogous formula to that for {W^{T)) 
in Eq.® but for the replacement of one of the operators 
bFf by Wj 4 yields the contribution of to the inter¬ 
action energy. Upon applying the operator — we 

find, at order 0{v/c), 


Q 2 /*-i-oo 

/T^fi-cx ^ —ZQ/C A B B A i In / 7 2 

) —_^9 A l^q ^ Im'^ / dk k 

* J —OO 
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where the superscript i? — c in (F^“‘^) denotes ’Rdntgen- 
conservative’, eum is the three-dimensional Levi-Civita 
tensor and the causal-adiabatic approximation has been 
considered with 77 —>• O’*". As explained in Ref, [l^ . the 
latter amounts to considering averages in time over a 
time-interval much greater than |A)(]j|. Mathematically, 
this is achieved by replacing each interaction Hamilto¬ 
nian W by We^* within the time integrals, with 77 —>■ O’*", 
and by extending the lower limits of integration to — 00 . 
In the above equation g{kR) is the dyadic Green’s func¬ 
tion of the magnetic field induced at a distance R by 
an electric dipole of frequency ck in free space. Their 
components are 
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After performing the frequency integrals in the complex 
plane we end up with 
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Straightforward evaluation of the Green’s functions at 
kARr reveals that, in the near field, ~ 

\AAB\kART/c <C 1, whereas in the far field the ratio goes 
like ~ \Aab\I ckA -C 1. From this we conclude that the 
conservative Rontgen force is always negligible. 

As for the non-conservative {R — nc) force, 

(Ffl-nc) ^ a^(^(0)|U+(r)d^ X B(R4)U^(T)|^(0)) 

(19) 

corresponds to the time-derivative of the Rontgen mo¬ 
mentum. It vanishes at zero velocity and an analogous 
calculation to the one for {P^AOiv)) be performed for 
V 7 ^ 0, i.e., the CM d.o.f. must be promoted to quantum 
variables. The only differences w.r.t. the calculation of 
replacements of an operator Wa with 
—X B(Ra) and of — iVRy with dr- This procedure 
leads to a formula similar to that in Eq. (TTOl) , 
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where the causal-adiabatic approximation has been ap¬ 
plied with ?7 —>■ 0+. Same as for Eg. lfTTl) . the v-dependent 
Green’s functions within the integrand can be split into 
Doppler terms and lag terms linear in v. The Doppler 
terms yield null contribution to the Rontgen force. The 
lag terms of the dyadic Green’s functions are, in addition 
to the second term of Eg. (ITT]) from G[fc|RT —v(T —1")|], 
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a term coming from C/[A:|Rt — v(T — t")|] which is also 
proportional to V}i{T — t")/RT, 
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respectively, with vj^ being the component of v perpen¬ 
dicular to Rt- The reason for the neglect of in 

the calculation of is that its contribution does 

not survive orientational average there, whereas it does so 
in (F^~"°). Finally, performing the time and frequency 
integrals of Eq. (1201) we are left with 
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xRedr f dt f dt'f dt"e’'(‘+*'+*"){ 

X [lm{G”"”(fc'RT)}Im{a^«[fc|RT - v(T - t")l]} 
+ Im{aP«(fc'i?T)}Ini{G™”[A:|RT - v(r - t")|]}] 


where the prime on the big square brackets denotes the 
derivative w.r.t. k of the functions therein. Note that, 
in contrast to its cancellation found in Ref.@ for an 
atom in an incoherent superposition of states, vanishing 
of (F^”"'^) does not hold for an atom in a metastable ex¬ 
cited state. On physical grounds, a non-vanishing value 
of the Rontgen momentum must be expected for a po¬ 
larizable atom A in motion and subject to the electric 
and magnetic fields induced by a second atom B at rest 
[^ . Less obvious is the fact that the Rontgen momen¬ 
tum varies in time. 

Straightforward evaluation of the Green’s fnnctions 
in the equation above reveals that, in the near field, 
(F^-"c)/(F«iW^^) ^ k\R\. <C 1 and hence it is negli¬ 
gible. On the contrary, in the retarded regime, after per¬ 
forming the orientational average over the random orien¬ 
tation of the electric dipoles, we find 
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(22) 

whose radial component equals 

Finally, as we did for the velocity-dependent vdW force 
we compute here the Rontgen force on the excited atom 
A as it flies with velocity v parallel to a thin plate (R) 
of ground state atoms of the kind R, at a distance d. 
Assuming that kAd^ 1, we obtain 
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and two more terms proportional to v±{T — t")/R r com- The dominant term in this equation comes from the in- 
ing from the expansion of the unit radial vector within tegration over the plane of the component of (F^””'”) 























which is parallel to . Interestingly, this term is an or¬ 
der kAd or uiaI^ab greater than the velocity-dependent 
vdW force in the retarded regime. 

V. CONCLUSIONS 

We have performed the time-dependent calculation of 
the velocity-dependent forces acting on an excited atomic 
dipole which moves at constant velocity with respect to 
ground state atoms of a different kind, in the quasireso¬ 
nant approximation at zero temperature. To this aim, we 
have used time-dependent quantum perturbation theory. 
Both the interaction with a single ground state atom and 
with a thin plate made of a random distribution of in¬ 
dependent atoms have been evaluated. We find that, at 
leading order, these interactions are linear in the velocity 
and no relativistic corrections enter at this order. 

In either case, the total force consists of a conserva¬ 
tive van der Waals component and a non-conservative 
Rontgen component. For their computation, the posi¬ 
tion and kinetic momentum vectors of the CM of the 
atoms have been promoted to quantum variables. We 
have shown that relativistic corrections to the van der 
Waals force enter at order (u/c)®. 

In the non-retarded regime the van der Waals force 
dominates, being always antiparallel to the velocity - 


hence can be referred to as quantum friction. On physical 
grounds, this force corresponds to the velocity-dependent 
recoil experienced by the excited atom during the pro¬ 
cesses of absorption and emission of the photons that 
it exchanges with the ground state atoms, which takes 
place at a rate |A^b|- This force is indeed proportional 
to vuja,bIcAab and hence to the emission rate of the ex¬ 
cited atom which is induced by the ground state atoms. 

On the contrary, in the retarded regime the velocity- 
dependent forces oscillate in space, the van der Waals 
and the Rontgen forces are of the same order in the atom- 
atom interaction, and the Rontgen component dominates 
in the atom-surface interaction. The latter corresponds 
to the time-variation of the Rontgen momentum, which is 
also mediated by the periodic exchange of quasi-resonant 
photons. 
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